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I. INTRODUCTION 

The extent to which the normal state of the high Tc su- 
perconductors is anomalous has long been debated. The 
most conclusive evidence for a break-down of Fermi liquid 
theory has appeared relatively recently via a (pseudo)gap 
in the fermionic spectrum with onset temperature T*. 
Recent Nernst-'^'^ and ac conductivity^ experiments have 
led to some of the most exotic indications for this pseudo- 
gap which, it is claimed, appear in the form of "vortices 
above T^" . More generally (but not universally) one as- 
sociates the complex of pseudogap phenomena with some 
form of precursor superconductivity. While normal state 
vortices are most closely associated with the well known 
phase fluctuation scenario^, a primary goal of this pa- 
per is to address these same anomalous transport data^'^ 
within the alternative pair fluctuation scheme^. In the 
process, we present a natural extension of time dependent 
Ginsburg-Landau (TDGL) theory and associated trans- 
port coefficients which addresses higher temperatures T* 
well outside the usual limited range of applicability near 
Tc- It necessarily follows that the bosons appear as quan- 
tum rather than classical fields. 

The pair fluctuations which we discuss here have as a 
natural antecedent the Gaussian fluctuations of TDGL. 
Indeed, when pseudogap effects are weak (T* « T^) 
TDGL theory has been found to provide a reasonably 
goo d^'^'^ representation of the non-critical fluctuation 
regime in the high temperature superconductors. In 
this way the presence of "pre-formed pairs" has been 
demonstrated, albeit, in the narrow temperature window 
near Tc- Two important features of the cuprates sug- 
gest that BCS theory is not applicable in that dominant 
fraction of the phase diagram associated with the pseu- 
dogap phase. The anomalously short coherence length 
^ in the cuprates combined with the observation of a 
non-vanishing excitation gap present at the onset of su- 
perconductivity has motivated a number of authors^ to 
contemplate stronger than BCS attraction so that the 
bosonic degrees of freedom appear at a temperature T* 
which may be significantly larger than Tc. Below Tc the 
counterpart of these pre-formed pairs appears as non- 



condensed bosons. The resulting superconducting state 
is mid-way between BCS and Bose Einstein condensa- 
tion (BEC)^iii*i^. One can, moreover, establishi^, that 
this BCS-BEC crossover approach is closely related to 
(Hartree-approximated) TDGL theory. 

Bosonic degrees of freedom are, in fact, stabilized 
by going beyond the usua6^ Gaussian approximation. 
When their interactions are included at the self consistent 
Hartree leveli^ii^ a fermionic pseudogap necessarily ap- 
pears; as a result of this depletion of the density of states, 
there are fewer fermions to which bosonic states can de- 
cay. At the Hartree level the presence of bosons and 
the existence of a gap in the fermionic spectrum are two 
sides of the same coin. These pair fluctuation approaches 
differ from other boson-fermion modelsi^ of high temper- 
ature superconductors because here the bosons naturally 
dissociate into fermions, and similarly fermion pairs may 
recombine. 

Hartree approaches to pair fluctuations also differ from 
the widely discussed phase fluctuation scenario which fo- 
cuses on fluctuations in the order parameter, and under- 
lying Mott physics. In the pair fluctuation approach the 
emphasis is on the observed small ^. Because it is a mean 
field theory, this scheme is not appropriate to the critical 
regime and critical exponents here belong to the Gaussian 
rather than 3d— XY class. While Coulombic effects are 
profound for order parameter fiuctuations in the ordered 
state, here they are less important, for precisely the same 
reasons that they are omitted in traditional TDGli. As 
in BCS theory. Coulomb interactions may also enter the 
binding and unbinding of pairs insofar as they renormal- 
ize or even, for non-s wave superconductors^", stabilize 
the attraction between fermions. 

While there are microscopic theories^ which serve to 
justify TDGL, the related diagrammatic formalism can 
become prohibitively complex and as a consequence, the 
bosonic contributions to transport properties are more 
readily deduced^ using TDGL theory directly. This, 
then, leads to the present focus on establishing a quan- 
tum extension of TDGL which is thereby amenable to 
detailed transport studies, at higher temperatures T*. 
To proceed, we construct a simplified model of charged 
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bosons subject to quantum dissipation, which has TDGL 
(with a non-BCS parameter set) as its special hmit. Pair 
fluctuations involve the continuous dissociation and re- 
combination of fermion pairs. To simulate this behavior 
we follow Caldeira and Leggett and treat the fcrmions as 
reservoir harmonic oscillatorsiStiS coupled to the bosons. 
Coulomb and hard core interactions between bosons are 
treated at the same level as in Hartree-approximated 
TDGL. 

We end by summarizing the key experiments we will 
address here. We focus on those experiments where 
the bosonic contributions dominate over their fermionic 
counterparts. It is convenient to define the electrical cur- 
rent J° = a'E+a{ — VT). In contrast to the behavior in a 
Fermi liquid, a sizeable Nernst signal ly (reflecting a com- 
bination of components of the a and a tensors?-) appears 
at an onset temperature appreciably above Tc, called T*; 
this temperature is loosely associated with T* . In the in- 
sulating phase, the thermoelectric response function axy, 
as well as i' vanish at T = 0. Precursor effects are simi- 
larly observed'^ in the imaginary component, (72, of ct at 
100 GHz, which set in below T*, but appreciably above 
Tc- Thus far, these onset temperatures are significantly 
below their counterparts for axy Moreover, for the range 
of frequencies measured, these complex conductivity data 
can be fit'^ to a rescaled Kosterlitz-Thouless form which 
allows extraction of a phase correlation time t'(T). In 
view of the above a and a experiments, it has seemed 
rather mysterious that the diamagnetic magnetization?,, 
has a relatively non-existent normal state precursor. 

Finally, the dc resistivity has been measured somewhat 
systematicalljiSLSSiSS in the underdoped regime. There is 
evidence for an onset temperature T* of roughly the same 
order of magnitude as in the other transport experiments. 
At this temperature an enhancement of the conductivity 
is evident. However, the magnitude of the deviation from 
the fermionic contribution is not noticeably large except 
in the immediate vicinity of Tc- Analysis of these data do, 
however, depend on assigning a particular temperature 
dependence to this background fermionic contribution. 



II. THE UNCORRELATED BOSON MODEL 

A. Revisiting Time Dependent Ginsburg Landau 
Theory 

The time dependent Ginsberg-Landau equation 
(TDGL) 



f-^ ( Ai( A HV-gA(M)r 

+ ao(r)^(x, t) + a'l^PiK, 01 V(x, t) + D{x, t) (1) 

is a natural starting point for characterizing the dynam- 
ics of bosons. This equation describes the fluctuational 
regime of conventional superconductors in the vicinity 



of Tc, where the bosonic degrees of freedom are repre- 
sented by the wavefunction i/)(x, t). Here \q\ = 2e rep- 
resents the bosonic charge, {ip, A) is the electromag- 
netic potential which determines the EM fields via the 
usual formulae E(x, t) = — V(/3(x, i) — ^A(x, i), and 
B(x, t) = V X A(x, i). The units adopted in this paper 
correspond to SI, with /i = eg = fcs = 1. The dynamics is 
importantly controlled by D{x,t) which is a white noise 
function satisfying 



{D*{x't')D{x, t)) = 2j2Td{t - t')d{x - 



(2) 



Generally we contemplate complex 7 with 72 = Im7 and 
71 = Re7. 

In this paper we focus on a description in which bosonic 
degrees of freedom may be treated as uncorrelated. For 
the TDGL case this corresponds to applying a Hartree 
approximation to Eq. (^) so that, 

{■^ ( a HV-gA(x,t))^ 

+ a{T)i;{^,t)+D{^,t), (3) 

Here a{T) represents the (absolute value of the) bosonic 
chemical potential which vanishes at the superconduct- 
ing transition temperature Tc, which is to be distin- 
guished from T™^ at which = 0. Frequently one ig- 
nores the quartic term altogether, as in strict Gaussian 
fluctuation theories-^. An interpretation as well as jus- 
tification for TDGL can be provided by microscopic T- 
matrix approaches, based on the well known Aslamazov- 
Larkin diagrams^. With the inclusion of Hartree effectsi^ 
the fluctuations have characteristic Gaussian exponents 
while Hartree self consistency contributions lead to a 
slightly modified T-matrix schemsiSi. 

At the more macroscopic level of Eq.(^ the fcrmions 
in a superconductor are irrelevant. Nevertheless, micro- 
scopic theory makes it clear that the bosonic degrees of 
freedom correspond to fermion pairs; moreover, it is the 
fermions which are ultimately responsible for the complex 
noise parameter 7. The self consistent Hartree approxi- 
mation introduces a fermionic exitation gap^'* (or "pseu- 
dogap" ) which is present at the onset of superconductiv- 
ity. This depletion in the density of states is responsible 
for the fact that the transition temperature (contained 
in a(T)) is lowered, relative to T™^. 

Despite its significant success in describing conven- 
tional superconducting fiuctuations, TDGL has known 
limitationsi^. In the context of understanding high 
Tc superconductors, one of the most serious of these 
is the necessity of introducing artificial cut-offs in the 
fluctuational spectrumii, often to depress the fluctua- 
tional contributions to transport. Other extensions of 
TDGL have been proposed which involve introducing 
a modiflcatioR-'' to the BCS-derived parameter set of 
Eq.(^, although for some experimental strict BCS the- 
ory appears to work quite well. 

What is most perplexing about high Tc superconduc- 
tors is the appearance of "pseudogap" effects with on- 
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set temperature T* . As T* progressively increases away 
from Tc, a strict BCS approach to fluctuation-based cal- 
culations of transport appears to be invalidatec^ . It is 
the premise of the present paper that the precursor super- 
conductivity of the pseudogap phase evolves continuously 
from the conventional fluctuation behavior seen, for ex- 
ample, in overdoped high Tc samples. While one does 
not expect TDGL to hold for T significantly larger than 
Tc, our goal here is to propose a natural extension of this 
theory appropriate for T well above T^. In this regime, 
the relaxation time of the bosons becomes comparable 
to or smaller than h/ksT: thus, the classical fluctuation- 
dissipation in TDGL should be replaced by a suitable 
quantum counterpart. In this way we treat tp in Eq.QJ 
as a field theoretic operator representing the annihilation 
of bosons. One of the most important parameters in this 
extended TDGL theory is T* which will enter into the 
bosonic chemical potential. This is the temperature at 
which the number of bosons vanishes. One may view 
T* alternatively as the onset of the fermionic pseudogap. 
These two viewpoints are two sides of the same coin, 
since bosons disappear or dissociate when fermions are 
no longer bound. 



B. Pre-Formed Pair Model: Extended TDGL 

To gain a deeper understanding of the essence of TDGL 
and of its prior success away from the pseudogap (or un- 
derdoped) regime^, we thus study a Hamiltonian describ- 
ing bosons on a d-dimensional lattice coupled to a quan- 
tum reservoir. Our treatment of the reservoir has strong 
similarities to the approach of Caldeira and Leggetlsi^. 
We consider 

H = ^euxV4(*) exp(-igrCux(i))V'x+u(i) 

UX 

X IX 



Here ipx{t) is the boson annihilation operator at lattice 
site X and time t, Annihilation operators for the reser- 
voir, Wi and Vi (with infinitesimal coupling constants r]i 
and Ci), are associated with positive and negative fre- 
quencies respectively, although the energies aj's and bi's 
are all positive. That two sets of reservoir operators are 
necessary will become clear later when we compare with 
standard TDGL. 

Here Eux is the hopping matrix element of the bosons, 
to be distinguished from its electronic counterpart, and 
Cux(i) = Jq dsUaAa{:!C+su,t). Surface effects appear via 
the X dependence of eox which also contains the Hartree 
interaction between bosons. 



The equations of motion of the system are given by 
— - q(p{x, t) j i/'x(t) = ^ £ux exp(-igCux(0)V'x+u(0 

' u 

+ J2vi^w,^{t) + 0x4x(*). (5a) 

i i 

d \ 

i— - qip{x,t)jvui^{t) = ai^Wi^(t) + ri*^ip^{t), (5b) 

'Wt " "-^^^ = -^-^'L(0 - QMt)- (5c) 



This model is manifestly gauge invariant. 

We solve Eqs. (|5bl) and l(5c|) in the temporal gauge 
(9?(x, t) = 0) to express Wi and v\ in terms of both i/'x(*) 
and their values at an initial time t^ — > —oo. Substituting 
into Eq. 15a|) . we find 

d 

ig^i'^t) = £ux exp (-igCux) '0X+U jt) 

U 

-if dt'j:2{t-t',T^)Mt')+D^{t), (6) 

where Tx is the local temperature and S2(i) the Fourier 
transform of 

E^(c^) = 2TrY,h\''Si^ - a^) - 2^^|C,P<5(c^ + b,), (7) 



It is reasonable to assume that this self energy, like any 
other, is smooth and that it vanishes as — s- ±oo. Here 

-Dx(i) = ^exp(-iaix(i - io))?7jxWix(io) 

i 

+ exp(+zfe,x(t - io))Gxwlx(io) (8) 

i 

is a function which represents a generalized or quantum 
noise. The coupling between the ip field and each reser- 
voir field is infinitesimal, so that the reservoir satisfies 
ideal Bose statistics and 

/I 
^S2(cc',rx)&(co,rx) 

•exp(-iw(t-t')), (9) 

with b{uj,T) = l/(exp(w/T) - l) the Bose function. The 
physics of the bosons, transport, magnetization, specific 
heat, density, etc, are governed by Eqs. l|njl and ©. 

The reservoir parameters a^, 6^, rji and Q are all sub- 
sumed into the boson self energy £2(0;). From this point 
forward we ignore these quantities in favor of the bo- 
son self energy. Moreover, £2(1^) depends only on lo 
for our localized reservoir. This simplification is sup- 
ported principally by the fact that this model captures 
the key physics found in microscopic scheme3^*2i, where 
it is sufficient to consider just the leading order k depen- 
dences. Differences between this previous strong coupling 
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T-matrix approach^^- and the present phenomenological 
boson model, are to be associated with the fact that the 
former scheme assumes that the T-matrix has fermionic 
constituents. These lead to a fine lo structure arising 
from the fermionic pseudogapSi, as well as to different 
high uj asymptotics. The ensuing simplification of boson 
physics, however, makes our model far more tractable. 

In the spatially uniform case, the Fourier transform of 
Eq. lO takes a simple form: 

(ij - £(k) + ^E;(co))v7(k,c^) = 5(k,c^), (10) 

with e(k) = Eu exp(iku), and 

%{u)^P I (11) 

We are now in a position to clarify the relationship 
between this theory and TDGL. When T is close to Tc, 
the bosonic relaxation rate is considerably smaller 
than T (typically of the order of tens or hundreds of 
Kelvin). In this regime the dynamics is dominated by 
low frequencies |aj| ^ T, and the Bosc function 

is well approximated by b{Lo,T) w T/lo. We presume 

that the high energy cut-off scale in I]2(w), which is 
associated with the reservoir, is of the order of typical 
electronic energies, thus thousands of K. Then, at low 
frequencies, the bosonic self energy is given by the linear 
functions 

e;(w)-i];(o)« (i-7i)w (12) 

and 

^2{i^) ~ 2j2CJ, (13) 

where from Eq. |(7J), 'S2{uj) vanishes at zero frequency. 

With the above approximations the quantum noise cor- 
relation function Eq. Q is reduced to its classical limit, 
given by Eq. . In the same way Eq. IjlUf) is reduced to 
Eq. with the important parameter 

-A*pair = a(T) = min£(k) + I];(0). (14) 

k 

We may now see that the two sets of reservoir fields, Wi 
and Vi, are associated with opposite charges, q and —q; 
they contribute to the positive and negative frequency 
regimes of S2(a;), respectively. If either of these two were 
omitted, S2('^) would vanish on one side of the origin; 
thus, the slope of S2('-^) would be discontinuous across 
zero frequency, and one would never arrive at the TDGL 
limit, no matter how small the frequency. In summary, 
the boson model presented in this paper is a natural quan- 
tum extension of TDGL. Conversely, TDGL is the low 
frequency limit of our quantum boson model. 

In preparation for computing transport coefficients, we 
introduce the bosonic correlation function, which can be 
derived from Eqs. H10() and (|5J. This is given by 

(?At(kV)?A(k,t^)) = A{k,u;)b{ij){2TTy+^5{k-k')S{u;-u;') 

(15) 



where A(k, uj) = Re2iT(k, uj) is the boson spectral func- 
tion, and 

f{k,uj)=(uj-e{k)-T^i{uj)+'-ir2iuj))' (16) 

is the boson propagator or "T-matrix" . The average 
number of bosons per lattice site is 

/Adl. J,, _ 
A(k,^)6(c.), (17) 

where v is the cell volume associated with each lattice 
site. 



C. Transport Coefficients 

The electric current J° and heat current follow from 
charge and energy conservation: 

1— n ( / Pi \ ^ ^ t 

J:(x,i) = ^|(^z--g^(x,t)j V'(x,i)| 

• ^TO''euxCxp(-igCux(^))V'x+u(i) -(-h.c. (18) 

u 

When a magnetic field Bah = gf~^& " applied 
to the system, surface electric and heat currents appear 
in a thin shell around its boundary. 

To derive these surface magnetizations we proceed as 
follows. We confine bosons within the sample boundaries 
via a spatially dependent eox which approaches the bulk 
value well away from the surface and -l-oo on the bound- 
ary. We presume that the spatial gradient of £ox is small, 
as is the applied magnetic field so that we may calculate 
the electric and heat currents to leading order in these 
quantities. Integrating the current in the normal direc- 
tion of the surface, from the boundary to deep within the 
sample, we obtain 

Ml, = i^i?,, y ^_A,,,,Rezr(k,c.)2c."6(c.), 

(19) 

where M^^ is the usual magnetization, is the ther- 
mal analogue^ of M^^, Aahcd(k) = ^ (wQc(k)w6d(k) 

-Uad(k)wfcc(k)), and i'afc(k) = g^^g^e(k) is the inverse 
mass tensor of the bosons, whose group velocity is given 
by va{k) = 4-e[k). 

When an external magnetic field is present, the bulk 
volume current J" must be combined with surface con- 
tributions so that the net "transport" currents are given 
bji^: 

Jl.^a = {J a) + -^El -f b,,.^Ml,El. (20) 

For small, but constant Bab-, electric field E'^ = 
Ea, and negative thermal gradient E^ = — g|-T, 
we obtain the linearized response functions J(l^)a = 
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Sri'=o X^b ^ab ' ^^'^ transport coefRcients 

L"" can be compactly written as 



T nn 
^ab 



^2 — n — n' 



2T"' 



6T«' 



d'^fc dw 
(27rp27r 

d'^fc dw 



^^^VaVcVMA'LU-+-'b('\ (21) 



(27r)'^ 27r 



where 6'^^^(a;) = — ^^j^ arises from the u> = limit, 
in which the boson absorbs or emits an infinitesimal 
amount of energy in the process of making a transi- 
tion to an adjacent energy level. = aat is the 
isothermal electric conductivity, L\l is the isoelectropo- 
tential thermal conductivity, and i"^ = aab and are 
off-diagonal coefficients. Eq. (|21|l satisfies the Onsager 
relation^ T""' L^]^' 0) = r"L;;^"(-B) (no summation), 
as a consequence of our inclusion of the surface termaSi. 
Surface effects enter into L]^ and Lj^, but cancel in 

rOO 
^ab- 

Finally, we deduce the boson contribution to the com- 
plex ac conductivity and ac Hall conductivity 



0'ab(ct;) 



.w„(k')«b(k')A(kV)&(t^') 



uj J {2itY 27r 

• (r(k', uo' + Lo) + r*(k', J ~uj)~ r(k'tj') - r*(k'tj')) 

■ (T(k', Lo' +oj)+f*{\i!,oj' -uj)- f{k'uj') - f* {k'uj')) . 

(22) 

It can be verified that Ijabi^ — > 0) = and that 
Eq. H22|) satisfies the f-sum rule. 



III. PHENOMENOLOGY 

In this section we arrive at a simple and generic phe- 
nomenology describing the pre-formed pairs or bosons 
of our extended TDGL model, in the context of hole- 
doped cuprates. Our goal in the next section is to use 
this phenomenology to address transport data at a semi- 
quantitative level with as few fitting parameters as pos- 
sible. It should be noted that this phenomenology is 
generally compatible with previously derived T-matrix 
based approaches^. The most important component of 
this phenomenological discussion lies in our introduction 
of the parameter T*, which has no natural counterpart 
in TDGL approaches. This is the temperature at which 
the number of bosons vanishes. This same temperature 
is reflected in the fermionic spectrum as that at which 
the excitation (pseudo)gap vanishes, so that fermions are 
no longer bound into bosons. The remaining parameter 
choices discussed below are reasonably straightforward, 
and represent natural extensions based on TDGL theory. 



There are three factors which govern the dynamics of 
our bosons: the dispersion function e(k), the self en- 
ergy S(w), and the chemical potential 

/^pair ■ 

We be- 
gin with £(k), which is given by £a(l — cos(fcaSa)) + 
^^(l — cos(fcbSb)) + £c{l — cos(fccSc)) + const., presum- 
ing only nearest neighbor hopping. Here Sc = s is the 
interlayer spacing, while Sa and Sh are in-plane lattice 
constants. The characteristic boson bandwidths {£a,b) 
satisfy £a ~ £b of the order of a few thousand K^^. By 
contrast, £c is on the order lOK for the least anisotropic 
cuprates (e.g., YBCO) and roughly two orders of mag- 
nitude smaller for the most anisotropic systems (e.g, 
Bi2212). In this way the dispersion function may be fur- 
ther simplified to yield fc^^/2TOa+i?c(l~cos(fccs))+const., 
the well known Lawrence Doniach dispersion. 

Next we turn to S(cj), which is expected to be given by 
its TDGL form. Eqs. (|12|l and for frequencies much 
lower than a characteristic cutoff energy fi, introduced 
earlier. Moreover, in the fermionic regime, far from the 
true bosonic hmit, (as appears to be the case for the 
cuprates), 71 <C 72- Our results are, thus, rather insen- 
sitive to 71 and for simplicity we set it to zeroi^ 

To obtain an estimate of O, we use our model to calcu- 
late the average number of bosons per lattice site n{T). 
At Tc we obtain the simple result'^^: 



n(Tc) « — ^ -f- — ^ In 



(23) 



where, for general energy scales E we define E = E/'-f2- 
It is reasonable to assume that well into the underdoped 
regime n{Tc) falls somewhere inside the range 0.01-0.5 
electron pairs per lattice site. Since Tc is considerably 
smaller than £a, the first term on the right side of Eq. H23|) 
dominates, and we deduce that is of the same order of 
magnitude as £a, thus thousands of K. The above anal- 
ysis indicates that the simple linearized expansion of the 
boson self energy, given by Eqs. (|12|l and H13|) . is a rea- 
sonably good approximation over the range of relevant 
frequencies we consider here. 

Finally we address the quantity /Zpair(T') which de- 
pends on the important temperature T* which has no 
natural counterpart in TDGL theory. This is the tem- 
perature where /ipair diverges. At this temperature the 
Bose degrees of freedom vanish. Concomitantly, at T* 
the fermionic excitation gap disappears, although, pre- 
sumably this is a crossover scale rather than a sharp 
transition. It has been argued that TDGL approaches 
overestimate fluctuation effects so that short wavelength 
cut-offs in the fluctuational spectrum have to be intro- 
duced, most recently to address the paraconductivitjiii. 
In the boson model introduced here, the non-TDGL tem- 
perature dependence in the pair chemical potential, /ipair, 
removes the necessity for introducing an ad hoc cut-off. 

We may estimate the magnitude of /ipair at intermedi- 
ate temperatures between Tc and T* , by using the boson 
density n(T). T-matrix based theories suggest^, that the 
electronic pseudogap scales with the number of bosons 
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above T^. A key assumption of our approach is that 
the magnitude of the pseudogap decreases by an appre- 
ciable fraction at intermediate temperatures between Tc 
and T*. This assumption is reasonably compatible with 
experiments such as ARPES, but it has not been con- 
clusively established at this time. We expect that the 
number of bosons evolves in a similar fashion. 

Within our theoretical model, we require that —//pair ^ 
in order to supress n{T) significantly relative to n{Tc), 
for T between and T*^. TDGL expressions for fipair 
(Mpair c>c (T — Tc)) wiU not lead to a sufficiently rapid de- 
cline in the boson number. It follows that there must be 
an additional term in the pair chemical potential which 
is negligible in the vicinity of Tc, but which rapidly in- 
creases at higher temperatures, and diverges at T*. A 
form generally compatible with the above analysis, as 
well as with microscopic T-matrix based schemes^, can 
be written as 



72 irri 



(24) 



where we have introduced a quantity 77 which represents 
the ratio of the boson coherence time t to its BCS value, 
7r/(8(T - Tc)), for T ~ Tc. Throughout, we take c - 
15, consistent with our order of magnitude estimates of 
Mpair- It should be stressed that there is nothing in this 
and the following sections which depends on the specific 
details of this functional form, provided the non-TDGL 
contribution (or second term) diverges at T* and vanishes 
sufficiently rapidly as T approaches Tc. 

Both ac and dc conductivity data indicate that 77 is 
close to 1 for optimally doped cuprates^^ but it becomes 
significantly larger than 1 as the material is progressively 
underdoped^. In this way the pairs live longer in the 
vicinity of Tc than expected from BCS theory. Quanti- 
tative analysis of the ac conductivity data (presented in 
the next section) shows that 77 can be of the order 10-20 
for a typical undcrdopcd Bi2212 sample. We take rj as 
the only essential fitted parameter in all of our numerical 
analysis^^ . Other parameter choices^ throughout this 
paper are for illustrative purposes only. 

We end this section by noting that there is an addi- 
tional quantum statistical effect which acts to further 
suppress boson transport. This is associated with the fact 
that the factor —b'{uj) is smaller than its TDGL counter- 
part T/lo^. This effect is appreciable when r ~ 1/T, 
and greatly amplified when r ^ 1/T, thereby further 
suppressing bosonic contributions to (DC) transport. 



IV. NUMERICS AND QUALITATIVE 
COMPARISON WITH EXPERIMENT 

A. Nernst Effect and Magnetization 

In Fig. n w6 plot the transverse thermoelectric co- 
efficient axy vs T for three underdoped samples with 
indicated T* /Tc. The arrows show the corresponding 



values of Tc which progressively decrease as T* /Tc in- 
creases. The inset plots (unpublished) experimental data 
on La2_j;Srj;Cu04 from Ref. which curves are for 
X = 0.12, 0.07, and 0.03 and which are in rough cor- 
respondence to the values of T* /Tc in the three theory 
plots. The onset temperatures are of the order of T*/2, 
as in experiment. In the vicinity of Tc, the calculated 
behavior of axy corresponds to that of TDGL theory, al- 
beit with modified coefficients. Thus axy diverges at Tc 
although V is finite there^. Similarly, with overdoping 
the behavior becomes characteristic of a more conven- 
tional fluctuation picture^: in which T* « Tc- The es- 
sential distinction between TDGL and the present case 
is that the onset temperature for bosonic contributions 
can be substantially higher than the conventional fluctu- 
ation regime. Semi-quantitative agreement between the- 
ory and experiment appears quite satisfactory. 

In both theory and experiment the dotted lines are for 
a non-superconductor with Tc = 0. To apply the present 
theory to the insulating case we presume that the flrst 
term on the right hand side of Eq. 124|) is a non- vanishing 
number which corresponds to the chemical potential of 
non-condensed electron pairs. This constant is chosen 
to be 24K in order to fit the maximum of axy{T) to 
its experimental counterpart. This figure illustrates the 
fact that axy (as well as v) vanishes at T = in non- 
supcrconducting samples. This derives from the behavior 
of the Bose function which approaches a step function of 
LD at very low temperatures. 



777^=7.4 
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FIG. 1: Theoretical curves for the normal state a^y with 
variable T* /Tc; experimental counterparts^^ in inset for 
La2-^Sr:^Cu04 at a; = 0.03, 0.07, and 0.12. Both the the- 
oretical and experimental values of a^y are normalized by 
0.02(B/T)V/(Knm), a number obtained earliest at 2Tc for 
an optimally doped sample. 



An interesting inference from the data™ is that, while, 
there is a considerable precursor effect for axy, the (or- 
bital) magnetization drops to its superconducting value 
only in the immediate vicinity of Bose condensation. 
This quantity [given as M^^ in Eq. (|19|) ] is plotted in 
the inset to Fig. |21 for the same parameter set as in 
Fig. ^ The fermionic background (measured experimen- 
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tally) is neglible, so that the bosonic contribution nec- 
essarily dominates. The sharpness of the Meissner onset 
(which clearly reflects Tc and not T*) can be attributed 
to the small ratio of the boson velocity to the speed of 
light (and the small size of the hyperfine constant). Sim- 
ilar results hold in TDGL-based calculations. 



B. AC Conductivity 



As in the data, here, Tg and Q! are deduced parameters. 
The inset plots t' = as a function of temperature. 
We find that Tg is roughly constant in T, in contrast to 
the data which finds this quantity to be decreasing as T 
increases. Nevertheless, the agreement with experiment 
for all three quantities plotted in Fig. O is quite good. 
In this way one might argue that key features of the ac 
conductivity data in Ref. S and attributed to KT physics, 
may equally well he explained by pre-formed pairs. 



In the main body of Fig. |21 are plotted the real and 
imaginary components of the ac conductivity as a func- 
tion of T for ix'/(27r) k, 100 GHz for two different values 
of T* /Tc, as in the previous figure. Both cti = Recr and 
(72 = Imcr are finite at Tc with a2 larger than tJi . Because 
of the relatively "high" frequencies (cj « 5K ^ fc)^, the 
associated frequency dependence is not the asymptoti- 
cally low u) limit, near Tc, where it would vary as l/^/bj. 
It follows that, just as in TDGL theory, dci/dT is finite 
while dCT2/dT diverges at Tc. We find that the magnitude 
of a is about twice the experimental value at T^. This 
prediction is the most notable difference between theory 
and experimenli^. 

The theoretically deduced onset temperature T*, for 
precursor effects in (J2 can be estimated by noting that 
the fermionic background is relatively negligible so that 
(72 becomes appreciable when it reaches a few percent 
of its value at Tc. This corresponds to onset tempera- 
tures which are factors of 2 or so closer to Tc than those 
estimated from otxy 



t> 20 




FIG. 2: Real and imaginary components of normal state cr 
(normalized by ctq = e^/s) at 100 GHz vs T. Solid and 
dashed curves have same T* /Tc as for counterparts in Fig.0 
Associated magnetizations are plotted in inset, with the ex- 
ternal magnetic field chosen as the unit. 

In Fig. 131 we replot the calculated ac conductivity (for 
T7Tc = 3.1), foUowinff the Kosterhtz-Thouless (KT) - 
based analysis of Refy- In the main figure we illustrate 
these KT fits. These are the basis of an interpretation 
of conductivity data in terms of vortices above T^. We 
present plots of = tan~^(cr2/o'i) and \S\ as a function 
of uj/^l' obtained from a{uj)/aQ = T^{T)S{uj /n')^'-^ . 




FIG. 3: Fits to Kosterlitz-Thouless (KT) scaling of the con- 
ductivity which can be compared to the analysis in Ref|3. See 
text for details. 



C. Resistivity 

In the main body of Fig. 0] we plot DC resistivity 
curves, p vs. T for an optimal and underdoped sys- 
tem. In addition to the bosonic contribution we have 
added a fairly generic^® linear-in-T contribution from 
the fermions, as the background term. The underdoped 
sample has the same parameters as in previous figures 
with T* /Tc = 3.1. For the optimally doped case we take 
T* /Tc ~ 1.45 and r] k, 1. It is clear that the transition 
region is wider for the underdoped system. The increased 
separation between T* and Tc, and the much elongated 
boson relaxation time, are responsible for this widening. 
Our optimal results compare favorably with those in the 
literature^^. For the underdoped case, the wider than 
experimentally observedSi transition regime that we find 
in p is closely connected to the behavior seen in the ac 
conductivity of Figure 2. And the latter is compatible 
with the experiments of Refj^- 

Dimensionality also plays an important role in deter- 
mining the behavior of the paraconductivity. In the up- 
per left inset of Fig.^ we plot the resistivity for a compa- 
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FIG. 4: DC resistivity p (normalized by pq = s/e^) vs. T, for 
underdoped and optimally doped samples. The underdoped 
sample has the same parameters as in previous figures. The 
optimally doped sample corresponds to r; = 1. In the insets 
are plotted the p and axy curves for an underdoped, more 3- 
dimensional, system like YBCO, with rj = 7.5, and T* jTc — 
3.1. The arrow in the lower right inset indicates Tc. 

rably underdoped, but more three dimensional material, 
such as YBCO. Increased 3-dimensionality (wo increased 
c-axis coherence length) clearly shrinks the appararent 
transition range of p. This can be traced to the reduced 
bosonic density of states in the low energy regime. To 
see how the change of dimensionality affects the Nernst 
coefficient, we present a plot of axy in the lower inset for 
this same sample. The onset for Nernst is roughly 40K 
above Tc- This contrasts sharply with that of p which 
is at most a few K above T^. Thus, the same sample, 
with exactly the same parameters, appears to have very 
different onset temperatures for resistivity and Nernst. It 
should be noted that if one looks closely at the plot of p a 
small downward curvature develops well above Tc, which 
can be traced to T* effects. 



V. CONCLUSIONS 

In this paper we have investigated the effects of pre- 
formed pairs or bosonic degrees of freedom on normal 
state transport properties. Our goal was to provide in- 
sight into transport anomalies in the high temperature 
superconductors. These anomalies are associated with 
non- Fermi liquid- like observations, many of which seem 
to be rather similar to what is found in the presence of 
superconducting fluctuations. By contrast to the stan- 
dard fluctuations of TDGL theory, however, these pre- 
formed pairs appear at relatively high temperatures T* 
compared to Tc- Throughout our discussion, it should 
be stressed that we view T* as a "crossover" tempera- 
ture, rather than a sharp phase transition. The origin of 
these preformed pairs is not specified. They may arise 
from the strong attractive interaction which drives the 
superconductivity. (Indeed, attractive interactions in a 



non s-wave channel^ are possible, despite the presence 
of strong Coulombic effects). However, one may enter- 
tain as well other scenarios for the mechanism of pair 
formation. 

In the process of addressing transport we have devised 
a new formalism for extending TDGL into the quantum 
regime, well away from the transition temperature where 
the bosons condense. Our approach, moreover, allows 
exact calculations of all transport properties and, in the 
immediate vicinity of Tc, our results are equivalent to 
those of standard TDGL. Comparison with experiment 
is quite satisfactory. It should be stressed that no par- 
ticular fitting to the data was done, but rather in this 
paper we have explored the generic features of the model. 
Moreover, we have addressed a wide variety of different 
experiments: Nernst, ac and dc conductivity and diamag- 
netic measurements. Aside from introducing a quantum 
statistical treatment of the bosons, a key feature of our 
approach is the introduction of the important tempera- 
ture T* . This is a new concept, not present in TDGL. 
Pairs now form far from the critical regime. But at the 
same time T* provides an ultimate cut-off. There are no 
bosons (or fermionic pseudogap) beyond this regime. 

Our paper shows that transport anomalies are com- 
patible with the presence of bosonic degrees of freedom. 
The same inference is made from thermodynamic and 
ARPES experiments which deduce the onset of a gap in 
the fermionic spectrum. Indeed, we view these broad 
classes of phenomena as two sides of the same coin. One 
of the key conclusions of this paper is that onset temper- 
atures for transport anomalies vary from one experiment 
to another. Moreover, all are considerably lower than T*. 
The boson contribution to transport depends sensitively 
on the pair life time t. Only sufficiently long lived bosons 
significantly contribute to transport. Hence bosonic con- 
tributions to transport are suppressed well before T* is 
reached. 

Ours should be viewed as an alternative to the vor- 
tex scenario or related phase fiuctuation picture for ad- 
dressing normal state transport anomalies. Indeed, it has 
been recently suggested^^ that phase fiuctuations alone 
cannot explain pseudogap phenomena. The reasonable 
agreement between the generic results of the present the- 
ory and experiment for all the figures provides support 
for a preformed-pair alternative to this vortex scenario. 
Non-condensed bosons are also present below Tc and, 
thereby, will lead to some continuity between transport 
coefficients across the transition. In extending this work 
to the ordered phase, however, it will be useful to find 
a relationship between bosons and vortices. Along these 
lines is the dual representation explored earlier— in which 
vortices are the basic particles, rather than Cooper pairs. 

Work supported by NSF-MRSEC Grant No.DMR- 
0213745. We thank Y. Wang and N. P. Ong for unpub- 
lished data, and I. Ussishkin and G. Mazenko for very 
useful conversations. 
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der term in 71 is important, since the zeroth order term 
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For the underdoped samples, illustrated in Figs. I1I2I and 
and the underdoped curve in the main body of Fig. |1| 
we choose ?7 = 20 to fit the onset temperature of Imcr at 
lOOGHz^ 

Here the quantity that is determined is actually Eajyi- A 
similar situation holds for Ec. and /ipair. Due to the approx- 
imate linearity of S(lj), only combinations such as /ipair/72 
or m72 appear in low frequency (^ f2) observables. 
This formula assumes that 72 < Ea and Ec <S T. For 
larger values of Q,, the first term still increases with the 
ratio Q./Ea, but slower than ^/{A-kEo)- For Ec comparable 
to or greater than T, the second term will decay slower 
with increasing Ec, and won't be zero or even negative, as 
this small Ec formula would suggest. 

The pair hfe time is thus on the order of 1 /Q. at {Tc+T*)/2. 
If l/Jl ^ l/T, as seems to be the case according to our 
analysis, this implies that bosonic transport is negligible 
well before T* is reached, as is confirmed by all transport 
experiments including the Nernst effect. 
We write Tc = max(1.6r* - O.ST*^ - 0.18,0) in units of 
roughly lOOK, to semiquantitatively fit the phase diagram 
of hole-doped cuprates. For all the 2D systems, we take 
Ec ~ 0.05K, while for the 3D system illustrated in the 
insets of Fig. |11 we take Ec ~ 15K. These are chosen to 
be consistent with the observed 3D-2D crossover tempera- 
tures of various materials. 



